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Strong mean field dynamos require supercritical helicity fluxes 
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Abstract. Several one and two dimensional mean field models are analyzed where the effects of current helicity fluxes 
and boundaries are included within the framework of the dynamical quenching model. In contrast to the case with periodic 
boundary conditions, the final saturation energy of the mean field decreases inversely proportional to the magnetic Reynolds 
number. If a nondimensional scaling factor in the current helicity flux exceeds a certain critical value, the dynamo can operate 
even without kinetic helicity, i.e. it is based only on shear and current helicity fluxes, as first suggested by Vishniac & Cho 
(2001, ApJ 550, 752). Only above this threshold is the current helicity flux also able to alleviate catastrophic quenching. The 
fact that certain turbulence simulations have now shown apparently non-resistively limited mean field saturation amplitudes 
may be suggestive of the current helicity flux having exceeded this critical value. Even below this critical value the field still 
reaches appreciable strength at the end of the kinematic phase, which is in qualitative agreement with dynamos in periodic 
domains. However, for large magnetic Reynolds numbers the field undergoes subsequent variations on a resistive time scale 
when, for long periods, the field can be extremely weak. 
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1. Introduction 

Astrophysically relevant dynamos tend to have boundaries 
or are at least confined. In all practically relevant cases they 
are certainly not homogeneous. Exceptions are dynamos on 
the computer where triply-periodic boundary conditions are 
used. Despite such dynamos being so unrealistic, they have 
played an enormously important role in revealing the nature 
of catastrophic a quenching (Brandenburg 2001, hereafter re- 
ferred to as B01, Blackman & Brandenburg 2002, hereafter 
referred to as BB02). This applies in particular to the case 
of helically forced turbulence. For reviews regarding recent 
developments see Brandenburg et al. (2002) and Branden- 
burg & Subramanian (2005a). However, some important con- 
clusions drawn from triply periodic simulations do not carry 
over to the case with boundaries. In this paper we discuss in 
particular the saturation field strength and focus on the mean- 
field description taking the evolution equation of current he- 
licity with the corresponding current helicity fluxes into ac- 
count. 

Homogeneous turbulent dynamos saturate in such a way 
that the total current helicity vanishes, i.e. (J-B) =0, where 
B is the magnetic field, J = VxB//i the current density, 
and fiQ is the vacuum permeability. (Here and below, (...) de- 



notes volume averages.) This is a direct result of magnetic 
helicity conservation in the absence of boundaries (B01). If 
the turbulence is driven at a scale smaller than the box size L, 
i.e. the forcing wavenumber fcf exceeds the box wavenumber 
k\ = 2ir/L (so fcf s> k\), it makes sense to use a two-scale 
approach. We therefore write B — B + b and J = J + j, 
where the overbar denotes an average field suitably defined 
over one or sometimes two periodic coordinate directions, 
and lower case characters denote the fluctuations. The con- 
dition of zero current helicity then translates to 

(J ■ B) = — (j ■ b) (no boundaries). (1) 

Together with the assumption that the large and small scale 
fields are nearly fully helical and that the sign of the helicity 
of the forcing is positive, we have (J ■ B) w —ki(B )///rj 
and (j • b) w fcf (b 2 )//io- The important conclusion from this 
is that, in the steady state, the amplitude of the mean field 
exceeds that of the small scale field, with 

fcf 



(B 



{b ) (no boundaries). 



(2) 



Moreover, for large enough magnetic Reynolds numbers the 
small scale magnetic energy is in rough equipartition with the 
turbulent kinetic energy, i.e. (6 2 }//xo « (pit 2 ), but see BB02 
for a more accurate estimate for intermediate values of the 
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2 DESCRIPTION OF THE MODEL 



magnetic Reynolds number. In the case considered in B01, 
where the scale separation ratio fcf/fci is either 5 or 30, the 
mean field energy exceeds the equipartition value by about 5 
or 30 - independent of the magnetic Reynolds number. 

The assumption of full homogeneity, which can only be 
realized with triply periodic boundary conditions, was an im- 
portant ingredient in arriving at super-equipartition fields. In 
this paper we discuss the more general case of non-periodic 
boundary conditions. We use here a mean-field approach to- 
gether with the dynamical quenching model (Kleeorin & 
Ruzmaikin 1982, Kleeorin et al. 1995), which proved suc- 
cessful in reproducing the homogeneous case as it was ob- 
tained using direct simulations (Field & Blackman 2002, 
BB02, Subramanian 2002). In the steady state without he- 
licity fluxes, the dynamical quenching model agrees with a 
modified catastrophic quenching formula which includes a 
term from the current helicity of the large scale field (Gruzi- 
nov & Diamond 1994, 1995, BB02). 



2. Description of the model 

In the mean field approach we solve the induction equation 
for the mean magnetic field B together with an evolution 
equation for the magnetic component of the a effect, 



— = V x (UxB + S 

at 



da 



M 



dt 



SB 



B 2 



C , 



(3) 



, (4) 



where current density is measured in units where fi = 1, rj is 
the microscopic magnetic diffusivity, i] t is the turbulent mag- 
netic diffusivity, i? m = rjt /rj is the magnetic Reynolds num- 
ber, £ is the mean electromotive force that includes, among 
other terms, a term proportional to a M B. In the following we 
adopt the numerical value k{/k\ — 5 for the scale separation 
ratio. The derivation of the cum equation is this normalization 
can be found in BB02 without helicity flux and in Branden- 
burg & Sandin (2004, hereafter BS04) with helicity flux. The 
connection with cvm = ^rj-b has been accomplished by 
using the definitions B 2 q = [i pu 2 ms and rj t = ^ru 2 ms , so 
r/(3^o/o) = B 2 q /r]t (see BB02). In fact, the evolution equa- 
tion for «m is therefore nothing else but the evolution equa- 
tion for the small scale current helicity, j ■ b. 

Indeed, j ■ b constitutes a possible source of an a effect, 
although it occurs usually in conjunction with the kinetic a 
effect that, in turn, is proportional to the negative kinetic he- 
licity, u> ■ u, where u> = V x u is the vorticity. The two 
effects together tend to diminish the residual a effect. We 
emphasize, however, that this does not need to be the case, 
and that even in the absence of kinetic helicity the magnetic 
a effect needs to be taken into account. One example is the 
case of a decaying helical large scale magnetic field, where 
j b is, being generated from the large scale field. The asso- 
ciated «m acts as to slow down the decay; see Yousef et al. 

(2003) for corresponding simulations, and Blackman & Field 

(2004) for related model predictions. 



The full expression for £ can be rather complex. For the 
present purpose we restrict ourselves to an expression of the 
form 



(a K + a M ) B + S x J rj t J 1 



(5) 



where q:k is the kinetic a effect, rj t is turbulent diffusion, and 
5 can represent both Radler's (1969) fi x J effect and the 
W x J or shear-current effect of Rogachevskii & Kleeorin 
(2003, 2004). (Here W = V x U is the vorticity of the mean 
flow.) The importance of shear, S (= dU y /dx in some of the 
first cases reported blow), and of the two turbulent dynamo 
effects, a and <5, is quantified in terms of the non-dimensional 
numbers 

Cs = J21 C a = J~, Cs = —. (6) 

VtM Vth r] t 

In the following we restrict ourselves to cases where either 
C a or Cs are different from zero. 

The general importance of helicity fluxes has been iden- 
tified by Blackman & Field (2000a,b) and Kleeorin et al. 
(2000). Here we use a generalized from of the current helicity 
flux of Vishniac & Cho (2001) flux, as derived by Subrama- 
nian & Brandenburg (2004), 



fSS 

a 



4>ijkBjBk- 



(7) 



Under the assumption that V • U = 0, we show in Ap- 
pendix A that 

4>ijk — Cyc tijlSlk, (8) 

where Sik — \{Ui,k + Ukd) is the mean rate of strain ten- 
sor and Cvc is a non-dimensional coefficient that is of order 
unity (see Appendix A). In the following we consider Cvc 
as a free parameter. It turns out that there is a critical value, 
C vc , above which there is runaway growth that can only be 
stopped by adding an extra quenching term. One possibility 
is to consider an algebraic quenching of the total a effect 
(a = % + Q!m)- Here we use a rough and qualitative approx- 
imation to the full expressions of Kleeorin & Rogachevskii 
(2002) by using 



o / 1 + g a B jB; 



(9) 



where we choose g a = 3 as a good approximation to the full 
expression (see BB02 for a discussion in similar context). For 
a completely independent and purely numerical verification 
of algebraic and non-i? m dependent quenching of a k and cum 
see Brandenburg & Subramanian (2005b). 

We expect the critical value C vc to decrease with increas- 
ing value of Cs- However, since Cyc should normally be 
fixed by physical considerations (which are uncertain), the 
possibility of a critical state translates to a critical value of 
Cs, above which "strong" (or "runaway") dynamo action is 
possible. This effect is similar to the W x J effect in that it 
only requires nonhelical turbulence and shear. However, we 
will show that, unless there is also current helicity flux above 
a certain threshold, the field generated by the W x J effect 
effect alone is weak when there are boundaries and when the 
magnetic Reynolds number is large. For the same reason, also 
a effect dynamos produce only week fields, unless the current 
helicity flux exceeds a certain threshold. 
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3 RESULTS 



Once the current helicity flux is supercritical, it is im- 
portant to make sure that om is spatially smooth. This is 
accomplished by adding a small diffusion term of the form 
K a V 2 Q:M to the right hand side of Eq. @. (Typical values 
considered below are n a = 0.02^.) 

We consider both one-dimensional and two-dimensional 
models. In both cases we allow for the possibility of shear. In 
the one-dimensional case (-L/2 < z < L/2) we allow for 
a linear shear flow of the form U = (0, Sx, 0), so the mean 
field dynamo equation is given by 



dB x 
dt 

dB„ 



0£y 

dz 



d 2 B,, 



dz 2 



— = bt> x + — — 
at dz 



v- 



d 2 B v 



dz 2 



and the current helicity flux is given by 

z ■ T c = \Cyc S{B X - B y ), 



(10) 



(11) 



(12) 



where z is the unit vector in the z direction. Note that, accord- 

2 2 

ing to this formula, assuming S > 0, and because B x < B y 
in such a shear flow, negative current helicity flows in the 
positive z direction. This is also the direction of the dynamo 
wave for «k > 0. [We recall that for ok > 0, negative 
current helicity must be lost to alleviate catastrophic quench- 
ing (Brandenburg et al. 2002).] We assume vacuum boundary 
conditions, 



B x = B y = (on z = ±L/2), 



(13) 



which implies that n • ^Fq = on the boundaries. This 
property can be regarded as an unfortunate shortcoming of 
the present model, because, although current helicity can ef- 
ficiently be transported to the vicinity of the boundary, it is 
actually unable to leave the domain. On the other hand, this 
type of boundary condition was also used in the simulations 
presented in Brandenburg (2005, hereafter referred to as B05) 
where the absence of closed boundaries clearly did allow for a 
significantly enhanced final field strength (see also Branden- 
burg et al. 2005). (It is also possible that in the simulations 
current helicity of the small scale field got lost because of 
numerical dissipation on the boundaries.) 

In order to allow for a finite current helicity flux on the 
boundaries we also compare with the case of an extrapolating 
boundary condition, 

' >! ~ (on z = ±L/2 fox i = x,y). (14) 



cz 



Thus, if e = we recover the standard vacuum boundary 
condition, Bi = 0. For e > 0, the slope of Bi(z) is such 
that Bi would vanish outside the domain on fiducial reference 
points, z = ±(1 + e)L/2. 

In the two-dimensional simulations we consider two dif- 
ferent cases. In the first case we use Eq. Jl 3I > in the verti- 
cal direction and periodic boundary conditions in the hori- 
zontal. In the second case we use perfectly conducting and 
pseudo-vacuum boundary conditions on the four boundaries 
in a meridional cross-section, just like in the correspond- 
ing direct simulations (BS04). The pseudo-vacuum boundary 
conditions are applied on what would correspond in the sun 
to the outer surface and the equatorial plane. 



In all cases the initial magnetic field is a random field 
of sufficiently small amplitude, so that the nonlinear solu- 
tions grow out of the linear one. We have not made a seri- 
ous attempt to search for solutions that only exist as finite 
amplitude solutions. Furthermore, we assume that initially 
«m = 0, which is sensible if one starts with weak initial 
fields. Again, we have not made a systematic search for finite 
amplitude solutions that might only be accessible with finite 
initial values of a M . 

3. Results 

3.1. Reference case with no shear 

We begin with the case of an a 2 dynamo where S = 0, so 
there is no shear and hence no helicity flux. The dynamo is 
excited for a > Tfrk\. A possible solution that is marginally 
excited and satisfies the boundary condition H3\ is given by 

(1 + cos kiz\ 
-smktz 1 ; (15) 

see Meinel & Brandenburg (1990) for the more general case 
of non-marginally excited (but still only kinematic) solutions. 

We have solved Eqs i ll Ob and ill 11 numerically using a 
third-order Runge-Kutta time stepping scheme and a sixth- 
order finite difference scheme. An example of a solution 
is shown in Fig. \l\ The results are displayed in Table ^ 
and Fig. |2 Both simulation data and mean field models are 
roughly compatible with the relation 

(B 2 )/B 2 cx i?" 1 (with boundaries), (16) 

that was first found analytically by Gruzinov & Diamond 
(1995) for the same boundary conditions Jl 31 . The simulation 
data shown in Fig. |2] supersede earlier results of Brandenburg 
& Dobler (2001) at lower resolution and smaller values of 

2 — 1/2 

R m where the scaling seemed compatible with B ~ i? m ' . 
However, in view of the new results this must now be re- 
garded as an artifact of insufficient dynamical range, so the 
correct scaling is given by Eq. dl6> . Furthermore, the simula- 
tion results give agreement with the mean field model if the 
dynamo number, C a , is somewhere between 3 and 10. This 
appears compatible with the fact that the scale separation ra- 
tio in the simulation is kf /k\ = 5, and that this ratio gives a 
good estimate of C a ; see BB02. 

Returning now to the description of the mean field calcu- 
lations, we note that at large values of R m the system shows 
relaxation oscillations where the sign of the field does not 
necessarily change (so the period of the field agrees with the 
period of the rms value). The frequency given in the table is 

LO = 27r/T p eriod- 

The fact that the saturation field strength decreases with 
increasing magnetic Reynolds number is bad news for astro- 
physical applications. However, this result is in agreement 
with simulations, lending thereby support to the applicabil- 
ity of mean field theory. 

The situation may change by allowing the field to extend 
also in one of the perpendicular directions (say the x direc- 
tion), because then the field may have its main variation in 
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3 RESULTS 



3.2 Solutions with uniform shear 




Fig. 3. Field lines superimposed on a color/gray scale representation of the normal field component for an a 2 dynamo in two 
dimensions in the xz plane with C a = 3 and three different aspect ratios (1,2, and 4). Note that there remain only two cells 
even for large aspect ratios. 




10 12 3 



s 




Fig. 1. Field structure for an a 2 dynamo (solid line B x /B eq , 
dashed line B y /B cq ) together with om, for C a — 3, and 
i? m = 10 2 . 



this new direction. In particular, if the field were to display 
Beltrami-like behavior in that direction, it might be more sim- 
ilar to the case of triply periodic boundary conditions. How- 
ever, as can be seen from Table|2j this does not seem to be the 
case. Even at large aspect ratios the wavelength of the field in 
the x direction remains of the order of the extent of the do- 
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Fig. 2. Dependence of (B ) / B 2 q on the magnetic Reynolds 
number for a run with open boundary conditions and no 
shear both for the mean field model and the direct simulation 
(squares connected by a line, with approximate error bars). 
The diamonds and crosses refer to mean field models where 
C a is 3 and 10, respectively. 



main in that direction (Fig. [5}, and the field amplitude still 
decreases inversely proportional to the magnetic Reynolds 
number (Table 0. 

3.2. Solutions with uniform shear 

For positive values of q:k and positive values of S there are 
dynamo waves traveling in the positive z direction. This is 



Table 1. Mean squared field strength of dynamically 
quenched a 2 dynamos. For time dependent solutions B is 
also averaged in time, and the frequency of the solution is 
given. Note in particular the inverse proportionality between 
(B ) and R m in the first three rows (marked by asterisks). 



Table 2. Same as in Tabled but for two-dimensional calcu- 
lations for different values of R m and different aspect ratios 

L x /L z . The inverse proportionality between (B ) and R m 
can best be seen from the last three rows (marked by aster- 
isks). The effect of changing the aspect ratio can be seen by 
inspecting the first two and the fourth row. 



Rm 


Cot 


(B 2 )/B 2 cq 


w/(7? T fc 2 ) 


10 1 


3 


1.35 x 10" 1 





10 2 


3 


1.35 x 10~ 2 





10 3 


3 


1.74 x 10~ 3 


0.28 


10 1 


10 


7.85 x 


0.85 


10 2 


10 


1.49 x 10" 1 


0.64 



Rm 

10 2 
10 2 
10 1 
10 2 
10 3 



3 
3 
3 
3 
3 



(B Z )/Bt 
1.06 x 10" 1 
2.30 x 10 _1 
1.25 x 10 _0 
2.53 x 10" 1 
2.71 x 10" 2 



1 

2 
4 

4 
4 
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3.2 Solutions with uniform shear 



3 RESULTS 



also the direction in which the flux of negative current helic- 
ity is pointing. Note that the saturation magnetic field strength 
decreases with increasing magnetic Reynolds number in the 
same way as before (Table[3j. The modified boundary condi- 
tion J14i with e 7^ does lead to an increase of the saturation 
field strength for R m = 10 2 , but not for R m = 10 3 . This 
behavior is not altered by the presence of helicity fluxes, i.e. 
changing the value of Cyc from to 0.2 has only a small 
effect. 

For Cs = 10, the critical value for runaway dynamo ac- 



tion is C 



vc 



0.15, but it decreases with increasing shear 



(e.g. for Cs — 20 we have C vc = 0.036.) This runaway 
growth may be a possible solution to the quenching problem 
in that it allows the solution to continue growing until it is 
saturated by other effects such as the usual a quenching that 
is independent of R m ; see Eq. 0. Note that when runaway 
occurs, the magnetic field increases sharply until it reaches a 
new saturation value close to equipartition; see Fig.|4] where 
Cyc = 1 has been chosen. However, the magnetic field is 
then no longer oscillatory. Obviously, the algebraic quench- 
ing adopted here is not realistic, but it does at least illustrate 
the point that when Cyc exceeds a certain critical value, there 
is runaway that could potentially be contained by having ad- 
ditional quenching terms. 

As was already anticipated by Vishniac & Cho (2001), 
a supercritical current helicity flux could by itself also drive 
a mean field dynamo. This mechanism requires a finite am- 
plitude initial magnetic field to get started; see Fig. [5] With 
unsuitable initial conditions the dynamo may therefore not 
get started and one might miss it. It is also possible that, even 
though the expected helicity flux is present and of the right 




0.00 0.02 0.04 0.06 0.08 0.10 
r) ki 2 t 



Fig. 4. Evolution of (B } / B 2 q for a model with Cyc = 1 
and additional algebraic quenching with g a = 3, n a = 
0.027/t, C s = 10, C a = 0.2, and R m = 10 4 . Note that the 
abscissa is scaled in resistive time units. 
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Table 3. Mean squared field strength of dynamically 
quenched a 2 fl dynamos. For time dependent solutions B 
is also averaged in time, and the frequency of the solution 
is given. For all calculations we used 64 meshpoints. For 
C a = 0.2 and Cs = 10 the kinematic growth rate is 

2 

0.08. The inverse proportionality between (B ) and R m can 
best be seen from the first three rows (marked by asterisks). 
Changing the sign of Cs has no effect on the saturation field 
strengths, regardless of the value of Cyc- 



Rm 


e 


Cvc 


(B 2 >/< 


Lo/(rj T kj) 


* 10 1 








1.52 x 10~ 2 


0.53 


* 10 2 








1.49 x 10" 3 


0.56 


* 10 3 








2.49 x 10" 4 


0.54 


10 1 





0.2 


3.18 x 10~ 2 


0.44 


If) 2 
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10 1 
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1.06 x 10~ 2 


0.56 


10 2 





-0.2 


1.35 x 10" 3 


0.57 


10 1 


0.2 





3.58 x 10" 3 


0.59 


If) 2 


0.2 





7.73 x 10~ 2 


0.36 


10 1 


0.2 


0.2 


4.31 x 10~ 3 


0.58 


If) 2 


0.2 


0.2 


1.45 x 10~ 2 


0.22 


10 3 


0.2 


0.2 


4.26 x 10" 4 


0.47 



Fig. 5. Evolution of (B ) / B^ q of dynamo with no kinetic 
helicity C a = 0, just shear with Cs — 10 and a supercriti- 
cal current helicity flux with Cyc = 1> an d different initial 
field strengths. (Because the initial field is random, much of 
the initial energy is quickly lost by dissipation in the high 
wavenumbers, which cannot be seen in the plot.) In all cases 
R m = 10 4 . Note that the abscissa is scaled in dynamical time 
units. 

kind, its strength remains subcritical (Arlt & Brandenburg 
2001). 

It is interesting to note that much of the late kinematic 
phase does not depend on the value of R m ; see Fig. [6] In 
fact, during the kinematic phase the runs with larger values 
of R m have slightly larger magnetic energies. Within a time 
scale that is independent of R m [here 2000 dynamical time 
units, (rjtkf the large scale magnetic energy reaches a 
significant field strength whose peak value is roughly inde- 
pendent of the magnetic Reynolds number. If one discards the 
subsequent decline of the magnetic energy, this result would 
be similar to that in the homogeneous case. In the weakly 
supercritical case, the cycle frequency does not decrease as 
nonlinearity becomes important. This is indeed an important 
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3 RESULTS 



3.3 Shear-current effect 




Fig. 6. Evolution of (B ) / £? 2 q for the a 2 SI dynamo with 
C a = 0.2, C's = 10, CVc = 0.2, and different values of 
R m . Note that the abscissa is scaled in dynamical time units. 



Fig. 7. Evolution of (B )/B^ q for the a 2 fl dynamo with 
C a = 0.2, Cs = 30, Cyc = 0, and different values of R m . 
The line styles are as in Fig. [6] above. 



feature of the dynamical quenching model that distinguishes 
it from the catastrophic quenching hypothesis (BB02). The 
run presented in Fig. [6] is for a finite current helicity flux 
(Cvc = 0-2), but the qualitative form of this plot is actu- 
ally independent of helicity fluxes and can also be obtained 
for Cyc = 0, for example. 

If the dynamo number is increased further to be highly su- 
percritical then the peak value reached by the large scale field, 
at the end of the late kinematic stage, increases even further; 
see Fig.0 This is again as anticipated by BB02 and Subrama- 
nian (2002). However, what was not anticipated are the sub- 
sequent strong dips in the mean field energy. It appears that 
the growth of ayi and the subsequent decrease of a, takes the 
dynamo below criticality. The mean field then decays until 
the microscopic diffusivity term [last term in Eq. @] causes 
cum to decay on a resistive timescale, and the net a to again 
increase, such that the dynamo becomes supercritical again. 
This qualitatively accounts for the long term oscillations of 
the mean field energy seen in Fig.0 whose period clearly in- 
creases with R m . It however implies that, on the average, the 
mean field energy again decreases with increasing R m , when 
there is no flux. 



3.3. Shear-current effect 

The issue of catastrophic quenching is often associated with 
the a effect alone, and it is implied that other large scale dy- 
namo effects may not have this problem. This is however not 
true. The main problem is quite generally associated with the 
helical nature of the large scale magnetic field. As is well 
known, the W X J effect can exist even without kinetic he- 
licity and just shear alone (Rogachevskii & Kleeorin 2003, 
2004). In the presence of closed boundaries, the current he- 
licity of the large scale field results in a corresponding contri- 
bution from the small scale field, which affects the resulting 
electromotive force. This can be seen quite generally by con- 




-3-2-10 1 2 3 








Fig. 8. Field structure for adxJ dynamo (solid line B x /B eq , 
dashed line B y /B cq ) together with for Cs = — 1, Cs — 
2, Cvc = 2, and R m = 10 3 . 



sidering the stationary limit of Eq. 10} for the case Cyc = 0, 
which yields a M — ~R nl £ ■ B / B^ q , 

So B 



SB 



(steady state). 



(17) 



1 + R m B IBl q 

Here we have defined the unquenched electromotive force 

£ a = a K B + SxJ-T] t J, (18) 

where the «m term is absent compared with Eq. (jSJi. To illus- 
trate the catastrophic quenching of the W x J effect we set 
a K = 0and<5 = (0, 0, <5) T . 

As is already clear from linear theory, a 8 x J effect can 
only produce self-excited solutions if S/S < (e.g. Bran- 
denburg & Subramanian 2005a), where the shear associated 
with W is responsible for stretching the poloidal field into 
toroidal. The 8 x J effect can also convert poloidal field into 
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3.4 Solar-like shear 



3 RESULTS 



Table 4. Mean squared field strength of dynamically 
quenched S 2 il dynamos. For all calculations we used 64 
meshpoints. The inverse proportionality between (B ) and 
R m can again be seen from the last two rows (marked by as- 
terisks). 



(-► 30° latitude) 



toroidal field, in addition to converting toroidal into poloidal, 
but this effect alone would not produce energy in the mean 
magnetic field. This is why shear is necessary. [We note, how- 
ever, that there is currently a controversy regarding the ex- 
pected sign of the S x J effect; see Riidiger & Kitchatinov 
(2005).] 

The field geometry is shown in Fig.|8j it resembles that of 
an a 2 dynamo (cf. Fig.Q, except that now both B x and B y 
are symmetric about the midplane, and cum is antisymmet- 
ric about z = 0. (The symmetry property of ckm is identical 
to that of £ ■ B, where one contribution is (8 x J) ■ B = 
V.B .) In Table 0] we present the results for different 
Again, we see quite unambiguously that, for 
a fixed value of S, the resulting field strength decreases with 
increasing R m , just like in all previous cases. 



values of R. 
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Fig. 9. Snapshot showing the field structure for an a 2 fl dy- 
namo with solar-like shear, C a = 3, Cs — 10 3 , Cyc = 
10 -3 , and R m — 10 2 . The solution remains highly time- 
dependent. 



3.4. Solar-like shear 

Finally, we consider the more complex flow geometry em- 
ployed by BS04 and B05 in an attempt to approximate the 
differential rotation profile in lower latitudes of the sun (cf. 
Figs 1 and 2 of BS04). In the simulations presented in these 
two papers the turbulence was forced in such a way that it has 
either positive, negative, or zero kinetic helicity. The latter 
case would correspond to no kinetic a effect, but the Vishniac 



Table 5. Mean squared field strength of dynamically 
quenched a 2 ft dynamos with a solar-like shear profile. All 

2 

solutions are time dependent, so B is also averaged in time, 
and the frequency of the field (not its energy) is given. The in- 
verse proportionality between (B ) and R m can be seen from 
the first three rows (marked by asterisks). For all calculations 
we used 128 2 meshpoints using C a = 3 and Cs = 10 3 . 
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& Cho mechanism and the W x J effect may still provide 
a possible explanation for the large scale field that is actually 
generated in such a simulation (B05). 

We use here the solar-like shear profile of BS04 that is 
given by 



U = (S*/fci)(0, cos kix cos kiz, 0) 



(19) 



in the domain — tt/2 < k\x < 0, < k\z < tt/2. As dis- 
cussed in BS04, k\x = —tt/2 corresponds to the bottom 
of the convection zone, where the toroidal flow is constant 
and approximately equal to that at 30 degrees latitude (corre- 
sponding to k\z = tt/2; the surfaces at z = and x — cor- 
respond to equator and outer surface, respectively, and permit 
current helicity fluxes. 

Length is measured in units of the inverse basic 
wavenumber of the domain, k^ 1 . Hereafter we assume k\ — 
1. The rate of strain matrix is then given by 

(0 sin x cos z \ 

sin a; cosz cos a; sinz , (20) 

cos a; sin z J 
so the divergence of the current helicity flux is, using Eqs 
and (|8jl, given by 



cos x sin z 



vc 



d_ l—i —2 

dx V 



B„ - B. 



SIM cos z 



7T ( B * B z 
oz v 



sin x sin z 



(21) 
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The corresponding value of V-^"p S is used in Eq. (0}, and the 
dynamo equation is solved subject to the same boundary 
conditions used in BS04. 

The field turns out to be highly irregular in time. An ex- 
ample of a snapshot at an arbitrarily chosen moment in time 
is given in Fig. |9] The magnetic field looks roughly similar to 
that found by averaging the results of direct simulations; see 
Fig. 6 of B05 or Fig. 18 of Brandenburg et al. (2005). 

In Table |5] we give the time and volume averaged val- 
ues of the squared mean field for different values of R m and 
different values of Cvc- Again, note that the energy of the 
mean field decreases inversely proportional to the magnetic 
Reynolds number. More surprisingly, increasing the value of 
Cvc has an adverse effect on the saturation field amplitude. 

4. Conclusions 

The present investigations still leave us with a puzzle. On the 
one hand the dynamically quenched mean field models yield 
invariably a resistively quenched saturation amplitude of the 
mean field - regardless of details of the boundary conditions, 
the presence of shear, or the nature of the dynamo effect. (In 
the absence of shear and just a 2 dynamo action, this result 
of mean field theory is well confirmed by turbulence simula- 
tions.) On the other hand, simulations with open boundary 
conditions (B05) have shown a clear difference compared 
with the case of closed boundaries. The reason for this dis- 
crepancy remains unclear at this point. It is however possible 
that Cvc is simply large enough, so there is a runaway dy- 
namo effect driving large scale fields by the Vishniac & Cho 
mechanism. A possible argument against this explanation, is 
that in B05, the dynamo did not resemble threshold behavior. 

Even though for subcritical helicity fluxes the saturation 
field strength may decrease with increasing R m , the field 
strength at the end of the kinematic regime seems to be still 
independent of R m . This is at least qualitatively similar to 
the behavior of homogeneous dynamos (BB02, Subramanian 
2002). Furthermore, the peak field strength at the end of the 
kinematic phase depends on the strength of the dynamo num- 
ber, which is also similar to what is predicted based on homo- 
geneous dynamo theory (cf. Subramanian 2002, Brandenburg 
& Subramanian 2005a). 

Clearly, further investigations of current helicity fluxes 
are warranted to pin down the origin of the apparently un- 
quenched saturation amplitude of the simulations. It was al- 
ready noted by BS04 that the Vishniac & Cho flux only ac- 
counted for about one quarter of the total current helicity flux 
that was determined from the simulations. This could also in- 
dicate that another perhaps more important component still 
needs to be included in the present mean field models. It is 
possible that the helicity fluxes discussed by Kleeorin et al. 
(2000, 2002, 2003a,b) may capture the missing components 
of the helicity flux, but this remains at this point only specu- 
lation. 

The present work has shown that there is a threshold 
of the current helicity flux above which runaway-type dy- 
namo action is possible. It remains a challenge to determine 



whether in fact all existing high magnetic Reynolds num- 
ber dynamos lie in this very same regime. Clearly, more 
work is needed to establish whether existing large scale dy- 
namos without kinetic helicity and just shear (B05) operate in 
this supercritical helicity flux regime by the Vishniac & Cho 
mechanism, or whether they work with the shear-current ef- 
fect, for example. 

Acknowledgements. We thank Eric G. Blackman for suggestions 
and comments on the manuscript. We also thank the organizers of 
the program "Magnetohydrodynamics of Stellar Interiors" at the 
Isaac Newton Institute in Cambridge (UK) for creating a stimulat- 
ing environment that led to the present work. The Danish Center 
for Scientific Computing is acknowledged for granting time on the 
Horseshoe cluster. 

References 

Arlt, R., Brandenburg, A.: 2001, A&A 380, 359 

Blackman, E.G., Field, G.B.: 2000, ApJ 534, 984 

Blackman, E.G., Field, G.B.: 2000, MNRAS 318, 724 

Blackman, E.G., Field, G.B.: 2004, PhPl 11, 3264 

Blackman, E.G., Brandenburg, A.: 2002, ApJ 579, 359 (BB02) 

Brandenburg, A.: 2001, ApJ 550, 824 (B01) 

Brandenburg, A.: 2005, ApJ 625, 539 (B05) 

Brandenburg, A., Dobler, W.: 2001, A&A 369, 329 

Brandenburg, A., Sandin, C: 2004, A&A 427, 13 (BS04) 

Brandenburg, A., Subramanian, K.: 2005a, PhR [arXiv: 

astro-ph/0405052 1 
Brandenburg, A., Subramanian, K.: 2005b, A&A [arXiv: 

astro-ph/0504222 1 
Brandenburg, A., Dobler, W., Subramanian, K.: 2002, AN 323, 99 
Brandenburg, A., Haugen, N.E.L., Kapyla, P.J., Sandin, C: 2005, 

AN 326, 174 
Field, G.B., Blackman, E.G.: 2002, ApJ 572, 685 
Gruzinov, A.V., Diamond, PH.: 1994, PhRvL 72, 1651 
Gruzinov, A.V., Diamond, PH.: 1995, PhPl 2, 1941 
Kleeorin, N.I. , Ruzmaikin, A.A.: 1982, Magnetohydrodynamics, 18, 

116 

Kleeorin, N., Rogachevskii, I., Ruzmaikin, A.: 1995, A&A 297, 159 
Kleeorin, N., Moss, D., Rogachevskii, I., Sokoloff, D.: 2000, A&A 
361, L5 

Kleeorin, N., Moss, D., Rogachevskii, I., Sokoloff, D.: 2002, A&A 
387, 453 

Kleeorin, N., Moss, D., Rogachevskii, I., Sokoloff, D.: 2003a, A&A 
400, 9 

Kleeorin, N., Kuzanyan, K, Moss, D., Rogachevskii, I., Sokoloff, 

D., Zhang, H.: 2003b, A&A 409, 1097 
Meinel, R., Brandenburg, A.: 1990, A&A 238, 369 
Radler, K.-H: 1969, Geod. Geophys. Veroff., Reihe II, 13, 131 
Roberts, PH., Soward, A.M.: 1975, AN 296, 49 
Rogachevskii, I., Kleeorin, N.: 2003, PhRvE 68, 036301 
Rogachevskii, I., Kleeorin, N.: 2004, PhRvE 70, 046310 
Rudiger, G, Kitchatinov, L. L.: 2005, PhRvE (submitted) 
Subramanian, K: 2002, Bull. Astr. Soc. India, 30, 715 
Subramanian, K., Brandenburg, A.: 2004, PhRvL 93, 205001 
Vishniac, E.T., Cho, J.: 2001, ApJ 550, 752 
Yousef, T.A., Brandenburg, A., Rudiger, G.: 2003, A&A 411, 321 

Appendix A: Vishniac- Cho flux in a shear flow 

We derive here the form of the Vishniac-Cho flux, introduced in 
Eq. Q, assuming that the underlying turbulence is homogeneous 
and isotropic (and weakly or not helical), and that all the anisotropy 
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required to get a non-zero <f>ijk arises from the influence of large 
scale velocity shear. We use the two-scale approach of Roberts & 
Soward (1975) whereby one assumes that the correlation tensor of 
fluctuating quantities (u and b) vary slowly on the system scale, say 
R. From Subramanian & Brandenburg (2004) we have 



Jspk 



where 



kikpV m s{k, R) d k, 



,(fc + ±K)u s {-k + \K).e lKR d 3 K. 



(Al) 



(A2) 



Here u m is the Fourier transform of the velocity component Mm. 
Note that for homogeneous isotropic turbulence 4>s P k vanishes. Now 
suppose we consider the effect of a weak shear on this turbulence. 
Then one can approximate its effect as giving an extra first order 
contribution to the velocity, u — it (0 ' + where u' ' is the 
isotropic homogeneous part. From the perturbed momentum equa- 



tion, u 



(i) 



r,[u m ■ VU + U ■ Vu {0) - Vp]. Here p is the 
perturbed pressure which ensures V ■ u' 1 ' = 0, and t* is some 
relaxation time. In Fourier space one can then write 



-,(!) 



(fc) 



T„P m j(fc) / [i44 0) (fc-fe')^(fe') 



+ i(k q - k' q )uf>(k - k')U q (k')] dV 



(A3) 



where P-mj^k) — $mj k m kj /fc 2 is the projection operator which 
ensures the incompressibility condition on u^ 1 '. Also U q is the 
Fourier transform of U q . 

We now substitute u = u (0) + u (1) in Eqs <A1> and IA2> . 
keeping only terms linear in u^ 1 '. Let us denote the two terms on the 



RHS of Eq. <A3j ; 



-(la) 



and u{ lb K Then on substituting Eq. <A3t 



into Eq. <AH . Eq. IA2> . four terms result, schematically of the form, 
Term I: Us > Term II: M m a ' ] u^\ Term III: u^ui lb \ Term IV: 
Urn^ Ms°' . The simplification of these terms involve tedious algebra. 
We outline the steps for Term I and simply quote the results for other 
terms. Term I is given by 



'->skp 



ATT*e klm / [kikpP m] (k + K/2)ik' q e lK - R UAk') 



X U, 



«>)/ 



{k+\K-k')u\ °'{-k + \K)\ d 3 Kd 3 kd 3 k'. (A4) 

We change variables to K' = K — k' and integrate over K', keep- 
ing in mind that the zeroth order u<°> is homogeneous. Also, since 
U only varies slowly with R, we retain only up to first derivative 
in U, which implies retaining only terms linear in fc' in the above 
integrals. One then has on evaluating the K' and fc' integrals 

<t>\k P = 4rr*efc! m V,{7 m / kik p v q °Jd 3 k 



4tt*A 
15~ 



^ksm^ p& m ~\~ €-klm V* I U jfi Sps 



(A5) 



Here we have taken the kinetic energy spectrum for the homoge- 



neous part of the turbulence to be v. 



P qa (k)E(k), and done 



the angular integrals over fc space. Also A — j E(k)k 2 d 3 k. Simi- 
larly we get for Term II: 



,H _ 4:TT„A | 

Vskp — I £fcs?7i V mU p + €klm V mU lOp. 



and so 



akp 



15 

+ 



/ii 

D skp 



= (8TT t A/l5) CskmSr, 



(A6) 

A similar calcu- 



I tv 

9 skp 



lation can be done for Terms III and IV to get 4> skp 
(8TTtA/3)e 3 kmSmp- Adding all the terms one gets the expression 
given in Eq. J8} of the main text, 



™ijk 



CvctijiSik; with CVc = 16tt*A/5. 



(A7) 



One can estimate the dimensionless number CVc as 

CVc = 16tt*A/5 ~ f (k c uT)(k c uT t ) ~ St 2 , (A8) 

where we approximated A = J E(k)k 2 d 3 k ~ ^u 2 k 2 , taken 
t ~ t» and defined a Strouhal number St = k B ur. For a flow dom- 
inated by a single scale k c ~ k{, the forcing scale. For a multi scale 
flow, like for Kolmogorov turbulence one should also keep the fc de- 
pendence of r (fc) oc fc -2 / 3 say, and fc e could be larger by a logarith- 
mic factor ln(fcd/fcf) where kd is the dissipative scale. However in 
a recent re-formulation of the dynamical quenching equation using 
local magnetic helicity density conservation (Subramanian & Bran- 
denburg 2005; in preparation) we have recovered the Vishniac-Cho 
flux as a magnetic helicity flux, and in this case fc c ~ k{. So it is 
reasonable to have St ~ ukfT < 1, and hence CVc < 8/5. 



SHeader : /home/brandenb/CVS/tex/mhd/helf luxld/pape 



1.98 2005/06/24 18:18:44 brandenb Exp S 



408 



